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Abstract 

We explore the use of bi-orthogonal basis for continuous wavelet 
transformations, thus relaxing the so-called admissibility condition on 
the analyzing wavelet. As an application, we determine the eigen- 
values and corresponding radial eigenfunctions of the Hamiltonian of 
relativistic Hydrogen-like atoms. 
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1 - Introduction 



Wavelet transforms have been successfully used in Mathematics, Physics and 
Engineering [Jl[ In particular, in the context of Quantum Mechanics, 

continuous wavelet transforms have proved very useful, giving rise - for ex- 
ample - to entirely new approaches to problems with spherical symmetry. 
For nonrelativistic Hydrogen-like atoms, an adequate choice of the analyzing 
wavelet reduces the radial Schrodinger equation to a first order differential 
equation, and the analyticity of wavelet coefficients leads, in a straightfor- 
ward manner, to the determination of the eigenvalues and their corresponding 
eigenfunctions ||. 

In this context, the selection of an analyzing wavelet is constrained by 
the "admissibility condition", which guarantees the existence of an inverse 
transform |5], |J. On the other hand, bi-orthogonal basis have been introduced 
in the context of discrete f| [7|] as well as continuous |§ transforms. 

However, in some cases, computational convenience may suggest that the 
most adequate "analyzing wavelet" be a non-admissible and even a non- 
square- integrable function. This is the case, for example, for relativistic 
Hydrogen-like atoms, as we will see later. 

It is the aim of this paper to extend the wavelet analysis to such situations, 
where it is not possible to construct an orthogonal continuous basis of L 2 , 
via the "ax + b" transform of the analyzing wavelet. In order to get an 
invertible transformation, we will rather restrict ourselves to a subspace of 
L 2 (containing the bounded eigenstates of the Hamiltonian to be treated) 
and make use of bi-orthogonal continuous basis. 

In section 2, we consider the space where this wavelet transform is well 
defined and some of its properties. We propose sufficient conditions for a 
function to belong to the space of wavelet coefficients. Such conditions are 
satisfied by the transformed eigenfunctions of the relativistic Hamiltonian to 
be treated later. 

In section 3, the radial Dirac equation is solved for relativistic Hydrogen- 
like atoms. As in the nonrelativistic case presented in reference ||, the 
analyticity of the space of coefficients is shown to determine the spectrum. 
Moreover, the aymptotic behaviour of functions in this space allows for a 
determination of the associated eigenfunctions. 

Finally, in section 4, we present some comments and conclusions. 
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2 - The transformation 

Let us consider a function ?/> 7 (g), solution of 

{lq + 2 ~J L ) r{q) = ~ r{q) ' (1) 
with q G [0, oo): 

ip-y(q) =g 7-2 e -g_ ( 2 ) 

For 7 > 1, ?/> 7 (g) is an admissible wavelet 0. So, by considering its 
"ax + 6" group transformation, 

^J(q) = a 3/2 e ibq [{aqf~ 2 e- aq ] , with z = b + ia and a > , (3) 

a continuous orthogonal basis of L 2 (R + , q 2 dq) can be defined as {ip]{q)} (For 
definetness, we will consider the radial part of three dimensional problems). 

Therefore, the wavelet coefficient of a function /(g) G L 2 (R + ,g 2 dg) is 
given by: 

(r z (q)J(q)) = ^- 1/2 F(z), (4) 

where 

F(z) = C< (/(g)) (z) = / dq e- iZ ^f(q) (5) 



is an analytic function of the variable z in the lower half-plane. One then 
has the reconstruction formula: 

o2 7 -2 . 

f(l) = T, pTo n\ I , <W*) /(?)) V>2(<2) (6) 

with dfiL(z) the left invariant measure of the "ax + 6" group: 
, . . dadb 

dfi L (z) = — — . (7) 

CL 

Moreover, the following equality holds: 

poo 2 7-2 „ 

/ dqq 2 \f(q)\ 2 = —— -/ dfJL L (z) {Im zY^\F{z)\ 2 , (8) 

Jo 27rr(27 — 2) J{imz>o} 
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which shows that F(z) belongs to a Bergman space 627-1 (see reference 

Now, for 1/2 < 7 < 1, the analyzing wavelet chosen is not an admissible 
one 0. So, in this range, it is not possible to construct an orthogonal basis 
leading to the reconstruction formula Moreover, for < 7 < 1/2, ^(q) 
is not even a square integrable function, so that the integral in equation (^) 
doesn't exist for an arbitrary f(q) G L 2 (R + , q 2 dq). 

In what follows, we will be interested in showing that it is still possible 
to use the transform in equation (|5]) for solving an eigenvalue problem, pro- 
vided certain regularity conditions are satisfied by its solutions. We will also 
analyze which properties of an authentic wavelet transform do still hold in 
such a situation. 

To this end, we will introduce a bi-orthogonal continuous basis. That is, 
we will make use of different functions in the process of analysis and later 
reconstruction: 

/(?)= / dfM L (z) M{q),f{q))ti{q)> ( 9 ) 

J{lmz>0} 

where xl(<?) is obtained - through the action of the group "ax + 6" - from a 
function x 7 (<?) satisfying: 

poo \ 

/ dqq^{q)*xi{q) = —. (10) 
Then, the following Lemmas hold: 

Lemma 1 Let f(q) E L} oc (R + , q^dq) H L 2 ((1, 00), dq)), with < 7 < I, and 
consider F(z) as defined in equation^). Then: 

a) F(z) is an analytic function in the half-plane lmz<0. Moreover, 
F(z) — > ; with Imz = a > 0, and F(z) —> 0. 

I He z\ — *oo Im z — ►00 

b) If f(q) ~ g a_1 (a > 0) for q ~ and f(q) is bounded when q — > 00, then 
C 1 transforms the operator qd/dq into the operator —zd/dz — (7 + 1). 

c ) If f(q) £ L 2 (R + ,q 2 dq) then: d^F(z) <G i3 2 ( 7 +i)-i, and 

POO 

/ dqq 2 \f(q)\ 2 = 
Jo 
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2 2( 7 +l)-2 



2vrr (2(7 + 1) - 2) Ji mz >o 



dfi L (z) (Jmz) 2(7+1) - 1 |^F(z)p 



(11) 



Proof: 

a) The function 

Fit) = F{b -ia)= / dqe- ibq q 1 f{q)e- aq 
Jo 

is the Fourier transform of q 1 f(q)e~ aq G L 1 (R + , dq). So: 



(12) 



F(z) 7^ , for Im z = —a < 

I Re Z\ — »oo 



(13) 



The analyticity of F(z) and the fact that F(z) — > for /mz — > —00 
are direct consequences of its definition (see equation(|12|)), since 
q 1 f{q)e~ aq G L^R+dg) for a > 0. 



b) Now, 

j\qe-^ q q 



/(<?) 



-iZg 7 



g 7 /(<?) 



(14) 



since the integrated term vanishes under the assumption made on the behav- 
ior of /(g), and e^ q q< +1 f{q) G L^R+dg). 



c) Notice that: 

d~ z F{z) = -1 / dqe- lZq qi +l f(q) = /7 +1 (/(g))(z) 
Jo 



(15) 



is the analytic factor of the wavelet coefficient of /(g) with respect to the 
wavelet ip] +1 {q) G L 2 (R + , q 2 dq), wich is admissible (since 7 + 1 > 1). Then, 
from equation (pi) we inmediately get equation (0). □ 
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Lemma 2 Let F(z) be an analytic function in the half-plane {Imz < 0}, 
with an asymptotic behaviour given by: 

F(z) =C (z- z )-^ +a) + C 1 (z- z- y^ +a+1) + G(z) , (16) 

where Imz > and \G(z)\ < K | ^| — (-r+ct+s) ^ \ oca Hy founded in the half- 
plane Imz < (C , C\ and K are constants). Then: 

a) {lmz)~ 1 ~ x l 2 F{z) is the wavelet coefficient of a function f(q) G 
L L(R + > (f'dq) n L 2 ((1, oo), dq), given by: 

/(?)=/ -/ db(Imz)^ 2 F(z)xi(q), (17) 
Jo a J-oo 

with z = b — ia. 

b) lfd~ z F{z) G B2 (7+ i)_i, then /(g) G L 2 (R+, g 2 ^). 

c) // |zc%Gr(z)| < X' |^|~( 7+Q!+2 ) ; and locally bounded in the half-plane 
Imz < (K' is a constant), then zd^F(z) = £ 7 (/i(g))(z) 7 where 

M?) = "(9^ + 7 + !)/(?)■ (18) 
Proof: 

a) In the first place, notice that (z — zo)~^ +a \ with a > 0, is the ana- 
lytic factor in the wavelet coefficient corresponding to the function fo(q) = 
C (> +a 7r( 7 + «)) q a - 1 e i *o q G L] oc (R + , g 7 dg) fl L 2 ((1, oo), dq). In fact, 

'Jo qq ~ [i(z - z- )p+" ■ U9j 

So: 

/ d^z) (/m^- 1 / 2 (z - ^o)- (7+a) Xl(g) = 

J{/mz>0} 

r oo , . jl+ a 

/ daa^ 1 x r {aq)2ir^- 1 [{z - ^r (7+a) ](<7) = w - , q a ~ l e^ . (20) 
jo 1 (7 + a) 
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(Notice that the integral in the first member is conditionally convergent). 
A similar result holds, changing a into a + 1, for the second term in equa- 
tion QISD, which is the analytic factor in the wavelet coefficient of fi{q) = 
d (> +a+1 )/r(7 + a + 1)) q a e^i G L[ oc (R+ q<dq) n L 2 ((1, oo), dq). 

As regards G(z), under the assumptions made, it belongs to the Bergman 
space -B27+1, since 

/ dfi L (z) (Im z) 2 ^ +1) - l \G(z)\ 2 < 00 , (21) 

J{lmz>0} 

as can be easily verified: For example, 

J J \G{z)\ 2 a 2 ^ 1 dbda < 00, (22) 

since |G(z)| is locally bounded. 

Moreover, is an orthogonal wavelet basis of L 2 (R + , q 2 dq), which 

defines a bijection onto B 2l+ i (see reference f|]). Then, g(q) G L 2 (R + ,g 2 <ig) 
exists such that: 

POO — 

G(z) = / dqq^e-^giq) = a +1 (g(q))(z) , (23) 
or, equivalently: 

fOO — 

G(z) = / dqqie-V'Mq) = F^fM^} = ^(f 2 (q))(z) , (24) 
Jo 

with f 2 (q) = qg(q) G L 2 (R+,<ig) and, therefore, f 2 (q) G L^R" 1 ", g 7 <ig) n 
L 2 ((l,oo),dg). 
Finally 

/ ^ L (^)(Jm^ 1 / 2 G(^)xl(g) = 

J{lmz>0} 
roo roo 

/ daaP , - 1 x r (aq) dbG(b - ia)e ibq = 

JO J -00 

2tt / daa^ 1 X J (aq)q^f 2 (q)e- aq = / 2 (g) , (25) 
Jo 

where use has been made of the fact that G(z) is the Fourier transform of a 
square integrable function. 
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b) Let us suppose that d z F(z) 6 -B27+1; then a function h(q) G L 2 (R + , g 2 <ig) 
exists such that: 



d z F(z) = / dqq 1+1 h(q)e~ iZq = C +1 (h(q)) (z) . 



(26) 



Moreover, from a), we know that: 



d z F(z) = d z / dqq>\f (q) + h(q) + f 2 (q)]e 



-iZg 



-iZq 



dqq" +1 [fo(q) + fi(q) + f2(q)}e 



(27) 



since the last integral is absolutely convergent. 

Then, from a) (with 7 — > 7 + 1), we conclude that fo(q) + fi(q) + /2(<?) 
f(q)^ih(q)eL 2 (R+,q 2 dq). 



c) In the first place, we will consider, for a > 0: 

/ dfx L (z) (Imzy- 1 / 2 \zd z (z - z o y^} X J(q) 

1 



— h+a) / dfi L (z) (Imz 

J{lmz>0} 



,7-1/2 



+ 



2f) 



(z - z )i+ a (z - z ) 7+a+1 



Xl(?).(2? 



From equation (|20|), the previous expresion reduces to: 
-(7 + «) 



? '7+ a - „-7+a+l 

ga-l e iZ Q q _|_ ^ ^(a+l)-l e i2:og 



r( 7 + a) 



T( 7 + a + l) 



,-7+a 



(g— + 7 + 1) 
ag 1 (7 + a 



(29) 



which proves the statement for the first two terms in equation fll6|) . 
As concerns the third one: 



dn L (z) (Imzy- 1 / 2 [zd z G(z)]x2(q) = 

db {b - ia)d b G{b - ia)e ibg 



{lmz>0} 

da a 7_1 x 7 (ag) 



roc r roo 

/ daa'- 1 x 1 {aq) - dbG(b - ia)(l + aq + ibq)e ibq 

JO l J-OD 



(30) 



where use has been made of the asymptotic behaviour of G(z) when inte- 
grating by parts. 

Notice that the integral between brackets in equation (|30|) is absolutely 
convergent, so that: 



dbG(b-ia)(l + aq + ibq)e ibq 

(l + aq + q-^j J°° dbG(b - ia)e lbq . 

Now, since G{z) G -B27+1, one has: 

1 r°° 

T~ x [G(b - ia)} (q) = — / db G(b - ia)t 

Z7T J-00 

with g(q) G L 2 (R + ,c/g). Therefore: 

+ aq + q-^j 2^g{q) e - aq = 

d 



ibq 



and 



2<Kq<e- aq U+ 1 + q-\g(q) 



di* L (z) {Imz)^ 1 ' 2 \zd- z G{z)]xl{q) 



{lmz>0} 

da a' y ~ 1 x y {aq) 



-2irqie- aq U + 1 + q j\g( q ) 



d 



(31) 



(32) 



(33) 



(34) 



which completes the proof. □ 



For < 7 < 1, the space of wavelet coefficients that appears in Lemma 
l.a) consists of functions (Jm2) 7_1//2 F(z), where F(z) is analytic for 
Imz < 0, vanishes for z — > oo and is such that d^F(z) belongs to £> 2 ( 7+ i)_i. 
This space of coefficients corresponds to the transforms of functions in 

LL(R + ,g 7 ^)nL 2 ((i,oo) )C /g). 

Now, we introduce the linear space A y of functions F(z), analytic in the 
half plane Imz < 0, vanishing for z — > oo and such that d^F(z) e £> 2 ( 7+ i)_i. 
Obviously, it is a pre-Hilbert space with respect to the scalar product: 

<F\G> A = [ dtiLiz^Imzf^+^dzFizYdzGiz) . (35) 

J lmz>0 



Lemma 3 The transformation C 1 , defined in equation (^) for < 7 < 1, 
maps a dense subspace of L 2 {H + ,q 2 dq) into a dense subspace of the pre- 
Hilbert space Ay, preserving the norm. 

Proof: 

Notice, in the first place, that the complete set of functions of L 2 (R + , q 2 dq) 
given by {ip n (q) = q a ~ 1+n e~ q ,n = 0,1,2...}, with < a < 1, is contained 
in L z 1 oc (R + , q^dq) H L 2 ((1, 00), dq). So, £ 7 is defined on a dense subspace of 
L 2 (R + , q 2 dq). Moreover: 

roo — 

C 1 (iffJq)) (z) = / dqq^ +a - l+n e- q e- lZq 
Jo 

= T( 1 + a + n) [i(z - i)]-^ +a+n ^ . (36) 
Now, the set {^(^(q)), n = 0, 1,2...} is complete in Ay, since 

id- z V (Mq)) (*) = / dqq^ +a+n e- q e- lZq = £ 7+1 (V„(g)) (z) , (37) 
Jo 

and because of the isometry established by the wavelet transformation £ 7+1 
between the Hilbert spaces L 2 (R + ,g 2 <ig) and £> 2 ( 7+ i)_i (see equation^)). 

Finally, for f(q), g{q) E L 2 (R+, q 2 dq) R (LL(R + , q'dq) R L 2 ((1, 00), dq) ) 
we have, 

<a(f(q))(z)\a(g(q))(z) > A = 
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<C' + \f{q)){z)\C^\g{q)){z) > B , 



27rr(2 7 - 2) 
2^~ 2 



(f,g) 



L 2 (R + ,q 2 dq) 



. □ 



(38) 



3 - Relativistic Hydrogen-like atom 



As an application of the results presented in the previous section, we pro- 
ceed, in what follows, to the determination of the bounded eigenstates of the 
Hamiltonian of relativistic Hydrogen-like atoms. 

As is well known ||11|| , after eliminating angular variables through the 
SU(2) symmetry enjoyed by the problem at hand, the radial part of the 
eigenfunctions satysfies the following equations: 



where m is the electron mass, and e are the allowed eigenvalues, satisfying 
|e| < m for bounded states. 

Moreover, A = Na (with N the number of protons in the nucleus and 
a = 1/137, the fine structure constant). In turn, \ is determined by the 
representation of SU(2) under study, and is given by: 





(39) 




(40) 



with j the total angular momentum of the electron. 
By defining: 



< q = 2r\/m 2 — e 2 , 
equation fl3"5|) can be rewritten as: 



(41) 
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te+H^-Gv^H™- - (42) 

where qf(q) and gg(g) are square- integrable. 

As it can be easely seen Hill , for q — > 0, the solutions of equation ([42]) 
behave as: 



/(g), <?(g)~<r 1+ ^% (43) 

with x 2 > A 2 . So, f(q),g(q) € L ; 1 oc (R + , q^dq), for 7 > 0. The transformation 
discussed in the previous section can therefore be applied since f(q) and g(q) 
satisfy the requirements of Lemma 1. 

Taking into account that the transformation is given by: 

F(z) = C? (f(q)) (?) = / dq e-^qy(q) , (44) 

J 

it is easy to see that (Lemma 1): 

£7 4=-(4 + H £7 ' (45) 

and: 

Oq = i^-C?. (46) 
dz 

So, transforming equations (S3), one gets: 



After some direct algebra, and calling 

•« = ( G« ) ■ < 48 » 

equation ([IT]) can be recast in the form: 

d I (A' + B' A'-B'\ _ . 

d=z^ z) =-2{^¥ + ^¥r zh m 
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with: 



M 



B' 



7 - X A 
-A 7 + X 



( A 



m+e 
m—e 



(50) 



(51) 



As is well known, the solution to equation (|49D is given by: 
1 



$(z) = Pexp<(-- /_ dz' 

2 Jzo 



A' + B' A'- B' 



z ' — - 

41 2 



(52) 



where P means ordering over the path leading from zq to z. 

Now, this expression can be greatly simplified through a judicious choice 
of 7: By taking [] 



7 = +VX 2 -A 2 >0, 
one has: 

(Mf = 2 1 M , MB' 



(53) 



2Ae 



,M 



m* — e c 



{B 



i\2 



2Xe 



y/m 2 — e [ 



,B' , B'M = 2~fB' 



and two new matrices can be defined as: 



.4 



A' + B' 
-277 



B 



where 



v = -7- 



A5 



i/m 2 — £ 2 



^4' — £?' 
-If) 



r] = -7 + 



As 



i/m 2 — s 2 



(54) 



(55) 



(56) 



1 Notice that, for x 2 = (j + 1/2) 2 < 1 + A 2 , 7 < 1, and we are in the conditions of the 
Lemmas of Section 2. 
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So, the following relations hold: 
A 2 = A , AB = A, 

B 2 = B , BA = B . (57) 
For this choice of 7 it is easy to see that ( |52"D reduces to: 
- $(z ) = 

Determination of the spectrum 

As discussed in Section 2, $(z) is an analytic function in the lower half-plane. 
So, its derivative: 



rf$ / rjA fjB \ (z 1 - |V (z' + i 



${z ) , (59) 



dz 2 ^ 2 / \ ^0 2 / V ^0 ~~~ 2 > 

must also be so. This requirement restricts fj to be a nonnegative integer: 

fj = —7 H — = n , n = 0, 1, ... (60) 

and 77 = —n — 27, from which the energy eigenvalues are seen to be: 

-1/2 

(61) 




A 2 



Vx 2 - A 2 + n 



Thus, as in the nonrelativistic case [f§], the bounded spectrum can be 
determined from the requirement of analyticity on the transform. 
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Determination of eigenf unctions 

From equation (p9|) and the condition $(z) — > for |z| — > oo (See Lemma 1 
of Section 2) it can be seen that: 

$(z) ~ z~ 2 \for \z\ -> oo. (62) 

So, the limit: 

lim ^ , = (63) 

is finite. 

Moreover, for 7 as given in equation fl53|) , the matrices A and S can be 
written as: 



m — £ T 



V m+e n / 



^Hral^)^ 1 ")' <65) 



Therefore, up to an overall multiplicative constant: 
$n{z) = (-7 + X- A i 



'm + e r 



m — e~ 



Z / „- x -(n+2 7 )-l , 



-7 + X + A 1 



m — e n \ + 



m-En 

m+e„ 



; / z "\ -(n+27) / ^ N n-1 



(OO 



dz'[z'--) + ■ (66) 



The integrals in equations (pq) can be evaluated on the imaginary neg- 
ative axis and analytically continued to the half-plane. In this way, one 
obtains: 



®n{z) = -7 + X - A 1 



m — e T 



m+e„ 



15 



;i + iz)-^ 2 Fx ( -n, 2 T , 2 7 + 1; — + 



-7 + X + 



'771 + E n 



m-e n J V +\/^+^ 



(1 + iz)~ 2 ~< 2 F, (-n + 1, 2 7 ; 2 7 + 1; — L^) , (67) 

where 2-^i(---) is a Gauss hypergeometric function. 

Notice that this solutions fulfill the hypothesis of Lemma 2, which guar- 
antees that we will obtain all the solutions in the configuration space. In 
order to do so, the "inverse" transform must be performed. To this end, an 
explicit x 7 (<?) function must be chosen. For convenience, we adopt: 

««) = SfH- (68) 

When inserting the first term in equation (|6"7| ) into equation fllTD , the integral 
to be solved is then given by: 

r N roc A n , rt 3/2 ibq , i 

lim db ^ a^ 2 (1 + iz)-^ 2 F X ( -n, 2 7 ; 2 7 + 1; 



n^ooJ-n Jo a 2 27rr(7) \ 1 + iz 

= — L-y(-l) fc H^2— i im [" dbe** daa?- l [l + ib + a]-^ +k) 

' e^q^e-i 1 F 1 (-n,2 1 + l;q), (69) 



T(2 7 ) 



where ii 7 i(...) is a degenerate hypergeometric function. 

The second term in equation ( |S7D can similarly be inverted (through the 
replacement —n — > —n + 1 in equation (|69"D). Thus, the eigenfunctions in 
the configuration space can be see to coincide with the well known result (as 
given, for instance in reference 
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4 - Conclusions 



In conclusion, we have explored the use of bi-orthogonal basis for continuous 
wavelet transformations, a generalization which is aimed at relaxing the so- 
called admissibility condition on the analyzing wavelet, and turns out to be 
useful for computational reasons. 

For definiteness, we have considered the radial dependence of functions in 
R 3 . As is well known, choosing as analyzing wavelet the function in equation 
(0), with 7 > 1, the wavelet transform in equation is an isometry between 
the Hilbert spaces L 2 {H + ,q 2 dq) and £> 27 _i. 

In Lemma 1, we have studied the transformation acting on functions 
f(q) G L} oc (R + , q 1 dq) fl L 2 ((1, oo), dq)), with < 7 < 1, a region where the 
analyzing wavelet is not admissible and can even be non square integrable. 
We have shown that the transform F(z) so defined is an analytic function 
in the half-plane Imz < 0, such that F(z) — >\R ez \^oo 0, with Imz = a > 0, 
and F(z) — >i mz ^oo 0, and that the transformation maps differential operators 
acting on f(q) into differential operators acting on F(z). Moreover, we have 
proved that, if f(q) G L 2 (R + , q 2 dq), then d^F(z) G £>2 7 +i- 

In Lemma 2, we have established that - for F(z) having an asymptotic 
behaviour as given by equation (|TB[) - the transformation has a right inverse 
through the use of a bi-orthogonal basis. 

In Lemma 3, we have shown that the transformation defined by equation 
(H), for < 7 < 1, is a mapping between a dense subspace of L 2 (R + , q 2 dq) 
and a dense subspace of a pre-Hilbert space Ay, which preserves the norm 
(defined in Ay in terms of the scalar product of derivatives in i3 2 y+i). 

Finally, as an example of the interest of our results, we have studied the 
spectrum of relativistic Hydrogen-like atoms. We have shown that, in the 
determination of eigenvalues of the Hamiltonian of this system and of their 
associated radial eigenfunctions, a wavelet transformation can be employed, 
and the calculation is greatly simplified by the choice 7 = +Vx 2 ~ F° r 
physical reasons, 7 can be any real number greater than zero, which makes ap- 
parent the need for our generalization of wavelet transforms. By applying the 
results proved in our three Lemmas, we have determined the spectrum from 
the requirement of analyticity on the transform, and we have reconstructed 
the associated radial eigenfunctions through the use of a bi-orthogonal ba- 
sis. Both the eigenvalues and eigenfunctions thus obtained can be seen to 
coincide with standard results. 
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